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TECHNICAL  REPORT  NO.  249 
IMPROVED  HYPERSONIC  LAMINAR  WAKE  CALCULATIONS 
INCLUDING  RATE  CHEMISTRY 

By  Martin  H.  Steiger 

I  INTRODUCTION 

This  report  is  concerned  with  the  prediction  of  the  properties  of  a  Laminar, 
axisymmetric,  hypersonic  wake  which  may  not  be  in  thermodynamic  equilibrium. 

The  flow  is  treated  as  a  continuum,  by  a  boundary  Layer  approximation,  utilizing 
integral  method  techniques  according  to  the  approach  developecyin  Reference  1. 

In  Reference  1  Bloom  and  Steiger, considered  both  laminar  and  turbulent  wakes. 
Their  integral -method  approach  involves  satisfaction  of  the  conservation  equations 
on  the  average  over  a  normal  section,  and  satisfaction  of  the  differential  equations 
exactly  along  the  central  axis.  Reasonable  profiles  in  terms  of  a  transformed 
normal  coordinate  and  (i  +  4)  undetermined  parameters  were  assumed,  where  i 
denotes  the  number  of  components  of  the  mixture  considered.  These  parameters 
are  utilized  to  express  the  axial  variation  of  the  physical  quantities  of  interest. 

The  (i  +  4  )  parameters  were  uQ  ,  HQ,  6m,  Qq  #  a^*  aNOo  anc*  the  e(luat^ons 
utilized  in  their  solution  were  the  momentum  integral,  momentum  boundary  condition, 
energy  integral  or  energy  boundary  condition,  and  conservation  of  species  i 
boundary  conditions.  The  development  of  the  solution  was  straight-forward,  but 
due  to  the  Limitations  of  a  small  digital  computer  used  for  the  computations  several 
approximations  in  the  calculation  procedure  were  made  which  may  have  caused  some 


error. 
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The  purpose  of  this  report  is  to  refine  the  aforementioned  work  in  preparation 
for  large -computer  caLcuLation.  In  particular  the  present  analysis  includes: 

(a)  An  improved  rate  chemistry  and  representative  thermodynamics. 

The  air  chemistry  includes  rate  processes  which  involve  seven  species,  i.e., 

Oz  ,  O,  N2  ,  N,  NO,  NO  and  e  .  The  pertinent  mass  action  laws  and  the  develop¬ 
ment  of  the  chemical  kinetics  and  the  thermodynamics  are  to  be  found  in  Reference  4. 

(b)  An  additional  parameter  is  introduced  (Hnn  ,)  which  couples  the  energy 
integral  equation  and  boundary  condition  •  This  permits  a  simple  evaluation  of  the 
effect  of  an  initial  defect  in  the  stagnation  enthalpy  and/or  Prandtl  -  Lewis  numbers 
different  from  unity  on  the  flow  field. 

(c)  The  system  of  equations  is  programmed  on  the  IBM  7090  digital 
computer.  This  enables  such  advantages  as  automatic  interval  change,  utilization 
of  a  Runge-Kutta  scheme  in  the  forward  integration,  test  for  maximum  truncation 
error,  evaluation  of  the  temperature  (which  is  given  in  terms  of  a  transcendental 
algebraic  expression)  to  any  desirable  accuracy  and  calculation  of  profiles  in  the 
physical  plane  at  any  predetermined  axial  stations. 

References  1  and  2  give  a  detailed  discussion  of  the  problems  associated  with 
wakes  behind  blunt  bodies  traveling  at  hypersonic  speeds  and  a  compilation  of  the 
more  recent  developments  and  references.  Of  particular  interest  at  the  present 
time  are  wakes  with  so-called  "turbulent  cores"  and  three-dimensional  effects. 

The  former  is  discussed  in  References  1  and  2  and  expanded  in  References  3  and  5, 
while  the  latter  is  considered  in  Reference  6. 

The  collaboration  of  Dr.  Martin  H.  Bloom  during  the  course  of  this  work  is 


acknowledged  with  thanks.  Also,  thanks  are  extended  to  Gertrude  Weilerstein  who 
programmed  the  analysis. 
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In  the  following  system  of  equations  for  mixtures  of  perfect  gases  it  is  assumed 
that  the  binary  diffusion  coefficients  of  several  species  are  approximately  the  same*, 
so  that  each  diffusion  velocity  Vj  is  given  by  Fickfs  law  =  -D  grad  a^.  There¬ 

fore,  a  single  Lewis  number  (Le  =  pDCp  /k)  appears,  wherein  the  fluid  properties 
are  those  of  the  mixture.  Likewise  a  single  Prandtl  number  (Pr  =  pCp/k)  of  the 
mixture  is  defined.  Both  parameters  subsequently  are  assumed  constant  along  the 
axis. 

The  equations  that  govern  the  flow  in  an  axially-symmetric  laminar  wake  are: 
Overall  Mass  Conservation 

(  pur  )x  +  (  p  v  r)r  =  0  (1) 

Momentum 

(  pu1  r)  +  (  puvr)  =  (  n  ru  )  (2a) 

x  r  * 

Pr  =  0  (2b) 


Energy 


(purH)x+  (pvrH)r  =  rHrJ  r  +  ^  ((T-l)ruur  j  r  + 


[flT  (Le-l)r  hi  <*ir] 


(3) 


*  For  an  elaboration  of  this  point  see  Reference  1,  Appendix  A. 


4 


Mass  Conservation  of  Species  i  (i  =  02  ,  N2  ,  O,  N,  NO,  NO  and  e  ): 
(purai)x  +  (pvrai)r=  r  airJ  r  +  P  r  Wi 


(4) 


where  notation  not  defined  in  the  text  is  defined  at  the  end  of  the  report. 

The  following  lists  the  boundary  conditions  (those  used  in  exact  solutions  are 
enclosed  in  curved  brackets)  and  auxiliary  conditions  employed  in  the  integral 
method.  The  latter  consist  of  the  conservation  equations  evaluated  along  the  axis, 
and  of  outer  edge  conditions,  which,  in  this  instance,  express  the  neglect  of  radial 
gradients  at  the  outer  edge. 

The  boundary  conditions  are: 

at  r  *J0  j  {  v  =0,  u_  =0,  Hr  =0,  =0  )  (5) 


:  {v  =0,  ur  =  0, 

puux  =  2  (  p,  ur) 


(6) 


P“Hx=2[f- 

Hr] 

r+  (<r-i)uur]  r  + 

+  2 

[f  ?  hi  air]  r 

r 

u  Le 

‘  air]  r  +  pWi 

pU  aix  =  2  1 

U  =  uQ(x), 

H  =Hq(x),  Hrr=HrrQ(x),  at 

U  =  ue,  H  ■  He,  a  i  =  Qie 

^where  ue, 

aie 

and  He  are  constantj 

ur  *  urr  “ 

o, 

H„  =  =0,  cu  =  aj  = 

r  rr  *vv 

(7) 

(8) 

(9) 

(10) 


(ID 
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Integral  equations  are  derived  from  Equations  (1-4)  assuming  that  integrated 
transport  terms,  that  is,  terms  like  ru^,  are  negligible  at  the  outer  edge  limit. 
The  integral  equations  are; 

Momentum 


p  u  (  1  -  u  )  r  d  r 


0=0  where 
c 


-  u)  r d  r 


(12) 


Energy 


or 

pu(l  -  H)  rdr 


(13) 


Species  Concentration 


_d 

dx 


pu  (  a i e 


-  Qt)  rd  r 


i_  ^  p  Wj  rdr 
ue  *4 


(14) 


The  following  transformation  are  now  introduced  into  the  integral  equations 
(  12  -  14  )  and  boundary  conditions  (6-8  ) 


m  d  m 

=  P  rdr,  m=6mn 

(15a) 

/6m 

62  /  2 

1  mdm 

=  J  — z: —  =  6 

0  n  m 

/  ndn 

Jo  7 

(15b) 

The  working  forms  of  the  governing  equations  and  boundary  conditions  are: 
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62  /  u(  1  -  u)  ndn  -  0 

m  -'o  c 


(16) 


uo  uo 


u 


X 


nn. 


e  m 


m 


(  1  -  H)  ndn  =  0, 


dH0  2ji, 


u. 


dx  p  6j? 

r  e  m 


Hnn0  (Oq.D 

+ 


C L 


CL 


+ 


(Leo  ”1) 


o  nn 


O  (L 


y;  hi  ai  1 

la  10  lnnQ 


(17) 


(18) 


(19) 


d 

dx 


[6m./0  «i>  "<*»].  ,6>  /  w,ndn 

m  o  i 


(20) 


d  a 


u 


i  o 


dx 


2Eq 

Pe6m 


ai  +  Wj 

'nn-  ' 


(21) 


and 


at  n  =  0,  u  =  uq(x)  ,  u  -■  0  ,  Hn  =  0 


(22a) 


‘i  , 


n 


H  s  H  (x) 
nn  nn-- 


n  -  1 ;  u  -  u 


H  -  *  a  .^  a  1 


(22  b) 


u  ~  u  ~  H  -  Ht,  ~  -  -  a  ,•  -  a ;  =-0 

n  nn  n  nn  ^nn 


A  compromise  between  simplicity  and  accuracy  influences  the  choice  of  specific 

procedures  for  solutions  by  the  integral  method.  Here  the  profiles  are  depicted  by 

polynomials  in  n  with  undetermined  parameters  (  u  .  H Hrn  and  a,  )  to 

o  u  nn0  a  q 

express  strearnwise  variation.  Another  undetermined  parameter,  the  wake  thickness 
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(6m)  is  included  implicitly.  The  equations  that  will  be  used  to  solve  for  these 
(i  +  4)  undetermined  parameters  will  be  Equations  ( 1 6  -  19)  and  (21).  Clearly, 
additional  parameters  can  be  included  and  additional  differential  or  integral  con¬ 
ditions  (for  example,  Equation  (20)  )  satisfied  at  will;  however,  the  more  usual 
approach,  underlying  simplicity,  will  be  followed  here.  The  profile  selection 
restricts  the  initial  profile  to  the  specified  form.  In  particular  cases  more 
accurate  curve -fits  of  the  initial  conditions  in  terms  of  n  can  be  employed.  Within 
reasonable  limits  these  will  not  have  a  substantial  influence  on  the  overall  results. 
The  peak  variations  represented  by  uQ,  H0  and  aiQ  are  believed  to  be  the  most 
influential  parameters. 

The  assumed  profiles,  which  satisfy  the  appropriate  boundary  conditions  are; 


=  6  n  ^  -  8  ri*  +  3  n  ^ 


(23a) 


and 


H  ®  HQ  +  (H0  -  Hq)  (10n3-15n4  +  6  n5 )  + 


2 


(23b) 
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III  CALCULATIONS 

By  utilizing  the  assumed  profiles  (23a-b)  the  governing  equations  (16-19) 
and  (21)  are  reduced  to  a  form  that  only  involves  the  (i  +  5)  undetermined  parameters 
and  variables  that  depend  explicitly  on  the  thermodynamic  state.  To  complete  the 
system,  these  equations  are  supplemented  by  the  thermodynamics  and  chemical 
kinetics  developed  in  Reference  4. 

This  system  has  been  programmed  on  the  IBM  7090  digital  computer  by 
Gertrude  Weilerstein* .  The  programming  techniques  in  the  present  case  are 
analogous  to  those  described  in  Reference  4.  Also,  the  present  program  incorporates 
a  procedure  by  which  both  the  normal  and  axial  distribution  of  the  pertinent  variables 
are  typed  out  at  predetermined  values  of  the  axial  coordinate.  The  calculation  pro¬ 
cedure  is  described  below. 

The  input  data  are, in  general,  the  conditions  at  the  edge  of  the  viscous  layer 
(i.e.,  ue,  Te,  pe,  Pe,  He  and  aj  ),  the  initial  data  at  the  axis  necessary  for  the 
forward  integration  (i.e.,  0C,  uoc»  H0£,,  aioc  an<^  ^nnQc)'  the  Lewis  and  Prandtl 
number  (LeQ,  (T0),  the  initial  value  of  the  streamwise  coordinate,  etc. 

At  the  initial  station,  say  x  =  xc»  sufficient  information  is  known  to  calculate 
the  following: 

(A)  Profiles  (for  0  <  n  <.  1  at  intervals  of  An  =0.1): 

u  =  u0  +  (ue  -  uQ )  (  6n^  -  8n  ^  +  3n4)  (24a) 

*  A  summary  of  the  data  preparation  for  the  7090  digital  program  has  been  prepared 
by  Gertrude  Weilerstein  of  GASL 
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^  4  c  ^nriQ  n  o  4  5 

H  =Hq  +  (He  -  Hq  )  (lOn  -15n  +  6n  )  +— j (n  -3n  +3n  -n  ) 

2  3  4 

ai  =  aio  +  (  a  i  e  -  ai  )  (  6n  -8n  +  3n  ) 

i  =  Oz  ,  N2  ,  O,  N,  NO,  NO+,  e" 

u2 

h  =  H  -  — 


(24b) 


(24c) 

(24d) 


and  for  the  corresponding  values  of  n,  the  temperature  (T)  is  calculated  from 


where 


and 


(24e) 


J  =  O 


2  • 


N 


2  » 


NO,  NO 


k  =  O,  N,  e 


i  =  Oz  ,  N2  ,  O,  N,  NO,  NO+,  e "  (24f) 


the  density  by: 
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the  physical  thickness  by: 

_[  210  ec  u2e 

(ue-uo)(10ue+lluo) 


Pe 

P 


(24g) 


and  the  particle  concentration  of  electrons  by: 

Qe-n  P 

Ne-=  — iJ-r  <24h) 

e 

Therefore,  profiles  of  the  thermodynamic  state  and  all  pertinent  variables  are 

calculated  as  functions  of  both  the  transformed  (n  )  and  physical  (r)  normal  coordinates. 

This  procedure  is  repeated  at  predetermined  values  of  the  axial  coordinate  (s). 

The  calculation  of  the  coefficient  cf  viscosity  at  the  axis  (  p.Q)  completes  the 

necessary  information  for  a  forward  integration.  It  is  given  by 

3/2 


p  =  5.48  x 


x  10 


-8 


(25) 


1.8  [T0+  199] 

The  axial  solution  is  then  obtained  from  the  following  system  (here  we 
define  s  =  x/L  where  L  is  a  characteristic  length  and  equals  unity  for  an  iaobaric 
flow  field,  therefore  L  =  1  in  this  analysis) 


d  u. 


L  Ho 


d  s 


Peue  35  uc 


(ue  -  uQ)  ( 10 u  +  11  u  ) 


(26) 


dH, 


ds 


4 

35 


LH, 


Pe0c  ue 


ue  -  uo 

u  „ 


(  10uo+  11  uQ) 


H 


nnQ 

TfcT 


(0-o-l)  (Lec-1)  y 

T  ^ 


uo (  ue  "  uo>  + 


i  hio*  °ie  ■°‘o) 


(27a) 
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where 


4*  hj  (a. 

i  Jo  i« 


-a 

lo  '  ie  IQ 


i0  )  =  ^hio  ai 


e  '  ho 


(27b) 


and  assuming  that  the  externaL  fLow  onLy  consists  of  the  undissociated  components  of 
air,  it  foLLows  that 

(27c) 


£hiQ  aie  =  hO,  aO-  +  hlvT  a 


O,  "Oz  -  “N2o  N2ft 


where 


and 


^2  0  7 

V  =  rt»  1 1  +  A°so 1 


N, 


hM  =  RT  _ iL  (I  +  A  .  ) 

N*o  °  Mm  2 


N2 


H  =  4(He-HQ)(349  ue  +  31 1  uQ)  _k 

nno  43  u  +  23  \i 

©  o 


4(-ue-uo)(10ue  +  11  u0) 
43  ue  +  23  uQ 


(28a) 


where 


K 


4(He  -H0c)(  349ue  +  311u0c)  -  Hnn0fl  (  43ue  +  23  u0<!) 


4(ue  -  u0c)(  10  ue  +  11  uCc) 


^2  8b) 


Obviously,  Equations  (27)  and  (28)  permit  an  estimate  as  to  the  effect  that  initiaL 
defects  in  stagnation  enthaLpy  and/or  Lewis -PrandtL  nurrjber  has  on  the  flow  fieLd. 
For  a  flow  which  is  initially  isoenergetic,  that  is,  Hfi  =  Hq^  and  HnnQ  =P,it  foLLows 
that  K  =0.  If  it  is  further  assumed  that  the  PrandtL  and  Lewis  numbers  are  unity 
then  Equations  (27)  and  (28)  state  that  the  flow  remains  isoenergetic  (i.e.,  H  =HQc= 
He  "constant  throughout  s  >  sc).  CLearLy,  under  these  conditions  the  governing 
partiaL  differential  equation  yieLds  the  same  results. 
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u 


dcti 


o  


4  ^eo 


ds 


Pe«*e«c  33  O’. 


*  ^(f  '■), 


(ue‘  uo)(10ue  +  Uu  )(  aj  -a,  )  + 

e  Ao 

(29) 


where  i  =  Oz  ,  N2  ,  O,  N,  NO,  NO  ,  e 

M4 

in  Equation  (29)  is  the  net  rate  of  production 


The  term 


w<)ic 


of  species  (  l/sec)  -with  contribution  from  i  mass  action  Laws  -  and  depends 
on  the  thermochemistry  of  air*  In  Reference  4  the  principles  c£  chemical  kinetics 
are  applied  to  the  following  reactions: 

02  +  X  - 


N,  +  X  ««- 


NO  +  X 


NO  +  O 


N  +  O 


NO  +  X 


20  +  X 

(30a) 

2N  +  X 

(30b) 

N  +  O  +  X 

(30  c) 

N  +  02 

(30d) 

NO  +  N 

(30e) 

NQ+  +  e' 

(30f ) 

NO+  +  e“  +  X 

(30k) 

(r  rjuuau 


Equations  (9a  w  k)  of  Reference  4  are  identical  to  the  expression 


for  i  =  Oz  ,  N2  ,  O,  N,  NO,  NO  ,  e  and  therefore  will  not 


be  reproduced  here. 
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Example  Calculation 

Calculations  have  been  made  of  the  laminar  wake  properties  for  one 
trajectory  point  of  a  blunt  body,  i.e.,  u^,  =  18,900  ft/sec  at  150,000  ft.  This 
corresponds  to  Case  1  of  Reference  1.  The  initial  data  is  shown  in  Table  1  and 
the  axial  distribution  of  the  pertinent  variables  in  Figure  la-d.  Several  typical 
type  -outs,  including  profiles,  are  given  in  Tables  (2a-f).  The  computation  time 
on  the  7090  digital  computer  was  approximately  five  minutes. 
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.c  specific  heat,  constant  pressure 

Dj  dissociation  energy  per  moLe,  (see  Reference  4) 

D^j  binary  diffusion  coefficient 

e  refer  to  electron  species 

h  static  enthalpy 

H  stagnation  enthalpy 

curvature  of  the  stagnation  enthalpy  profile  at  n  =0 
Ij  ionization  energy  per  mole  (see  Reference  4) 

k  thermal  conductivity 

K  defined  by  Equation  (28) 

L  characteristic  length,  this  report  L  =  1  ft. 

pDi_  c 

Le  Lewis  number  Le  = - *J  £ 

k 

m  transformed  normal  coordinate  Equation  (15a,.  b) 

M^  molecular  weight,  mass/mole 

n  transformed  normal  coordinate,  , Equation  () 5  a,  b) 

N  refers  to  atomic  nitrogen 

N2  refers  to  molecular  nitrogen 

Ne-  refers  to  particle  concentration  of  electrons,  electrons/cc 

O  refers  to  atomic  oxygen 

O*  refers  to  molecular  oxygen 

NO  refers  to  nitric  oxide 

NO+  refers  to  ionized  nitric  oxide 

r  normal  coordinate 

R  gas  constant  per  mole 


SymboLs  (Continued  ) 


s  transformed  streamwise  coordinate  s  =  x/L, 

T  temperature,  °K 

V 

characteristic  vibrational  temperature 
u,  v  streamwise  and  normal  velocity  component 
Uqq  free  stream  flight  velocity 


Wi 

x 

X 


a 


i 


6 


6 


m 


0 


0 


0. 

l 


P 

A 

S 

n 

0" 


net  rate  of  production  of  species  i 

streamwise  coordinate 

catalyst 

mass  fraction  of  species  i 
physical  viscous  layer  thickness 
transformed  viscous  layer  thickness 


momentum  thickness 


0 


“/pu(l 


energy  thickness  0^,  “  p  u  (1  -  H) 
,th 


itn  species  thickness  9  =  f  p  u( 

vo 

coefficient  of  viscosity 
density 

see  Equation  (24c) 
summation 


Le 


Avagadro's  number 
Prandtl  number  CT  ■ 


p  c 


-  H) rdr 
r  d  r 

-  ai)  rdr 
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c 


e 

e 


N 

N2 

NO 

+ 

NO 

O 

°2 

i,  J.k 


o 

x,  y,  n 
oo 


conditions  at  initial  station 

conditions  at  edge  of  viscous  Layer,  constant 

refers  to  electron  species 

refers  to  atomic  nitrogen 

refers  to  molecuLar  nitrogen 

refers  to  nitric  oxide 

refers  to  ionized  nitric  oxide 

refers  to  atomic  oxygen 

refers  to  molecuLar  oxygen 

indices,  (  i  =  #  ^2  *  ^0  >  NO  >  e  ^ 

values  evaluated  along  axis 

denotes  partial  differentiation  with  respect  to  indicated  variable 
undisturbed  flight  conditions 


Superscripts 

denotes  nondirnensiona L  quantities  with  respect  to  condition 

—  u 

at  the  edge  of  the  viscous  layer,  u  =  —  , 


) 


H 


etc. 
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A.  Flight  Conditions 

(1)  Velocity 

(2)  Pressure 

(3)  Density 

(4)  Temperature 

B.  Initial  Conditions 

(1)  Mass  Fractions 


ao 

- 

.  232 

ao2 

at 

0 

aN 

s 

.  133 

QN2 

= 

.  635 

(2)  Velocity 

(3)  9c 

<4)  Hhn0 

(5)  Stagnation  enthalpy 

C.  Miscellaneous 

(1)  Crocco  Integral  H 

(2)  Le0=  <T0  *  1 


1 8,  900  ft/ sec 
2.94  #/sq.  ft 

L  3 

3.39x10"  slug/ft 
280  °K 


a  NO  -  U 

a  ^,0+  =  5.63x10 

a  =  1. 029  x  10'8 

13,  078.  8  ft/sec 
.0065  ft1 2 
0 

1.  81646  x  10^  ft^/sec2 

constant  throughout 
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Figure  la.  Distribution  of  Temperature,  Velocity  and  Physical  Boundary  Layer  Thickness  Along  the  Axis. 
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igure  lb.  Distribution  of  and  q»,  Along  the  Axis 
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Figure  1  c.  Distribution  of  Enthalpy  and  Density  Along  the  Axis 
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Figure  Id.  Particle  Concentration  of  Electrons  A 1  on^  the  Axis 


